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Abstract 

A simple superfluid model with an effective four body interaction of monopole 
pairing type is used to explain the staggering of the charge radii in the isotope 
chains. The contribution of deformation and of the particle number projection 
are analyzed for the Sn isotopes. Good results are obtained for the staggering 
parameters and neutron pairing energies. 
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Detailed studies of the nuclear charge radii exhibit significant deviations from an average 
dependence of the mass number [I]. The systematic study of these deviation jl| revealed 
two important effects (1) the change in slope of the variation of the charge radii versus the 
neutron number when crossing a magic number and (2) the staggering of the charge radii 
in isotope chains (odd - even effect). These effects are found over the whole table of nuclei 
and are confirmed with increasing accuracy H - |J. 

The odd-even effect exhibit some important characteristics [§] (a) at the beginning of a 
neutron shell the amplitude of the effect is rising and it decrease to the end of the sell; (b) 
such variations in the amplitude are almost insensitive to the proton configuration; (c) at 
the end of a neutron shell an increasing of the amplitude for lighter elements is observed. 

Many attempts to explain the odd-even effect have been done during the last 25 years but 



the output was often contradictory. Uher and Sorensen |T(J have extensively calculated the 
charge radii for many isotopes chains, in the framework of pairing plus quadrupole model, 
allowing for monopole and quadrupole core polarization but, they were not able to describe 



the odd even staggering. Reehal and Sorensen ]TJJ were able to find agreement with the 
data for some cases taking into account the influence of the neutron blocking on the ground 
state quadrupole vibrations in a mixed, microscopic and liquid drop model. 

Calcium isotopes have been extensively studied ||, fll2| , the main mechanism involved 



being the core polarization, either in the Isospin Projected Hartree Fock (IPHF) framework 
or using a simple parameterization of the neutron proton interaction in a j n configuration. 
The IPHF mechanism [|L3[] gives only a qualitative agreement with the data but, the simple 
parameterization of Talmi |12| nicely describe the charge radii of calcium isotopes. However, 



the fitted parameters entering Talmi's formula have not an additional physical meaning and 
the realistic calculation of the core polarization contribution give a much smaller effect [0. 
The formal extension of the j n formula to the case of Pb isotopes has also received some 
criticism 0. 

Some peculiar cases, connected with large deformation effects, received particular atten- 
tion. The large amplitude of the staggering in neutron deficient Hg and Au isotopes has 



been attributed to an oblate - prolate instability |14|]. For the transitional nuclei Nd, Sm 



and Gd, it is possible to obtain, in some cases, a staggering of correct magnitude by taking 
into account the zero point motion of the nuclear surface and the effect of the h n / 2 neutron 



intruder state |0|. Other sophisticated calculations like the large shell model configuration 
for light nuclei [|nj or the HFB model taking into account the continuum states [T7| for Sn 
isotopes, have used to describe the odd-even effect. 

All the above mentioned models have not been able to put in evidence an underlying 
mechanism for this effect, observed with few exceptions on the whole table of atomic nuclei. 
If successful, they pointed out some specific contributions like the j n configuration, the 
prolate oblate instability, the effect of zero point surface vibration, etc. The decrease of the 
staggering toward the closure of the neutron shell indicate that the primary mechanism which 
discriminate between even and odd neutron number is pairing. This observation made by 



Zawischa |18j has been used to find a mechanism that strongly connect the neutron pairing 
properties with the proton properties. He looked for an universal mechanism which involves 
the collective properties of the neutrons and protons. The failure of the standard nuclear 
models to describe the staggering indicates that some small residual interaction can strongly 
contribute due to the collectivity. For the neutron system the collective quantity must be 
the pairing tensor, which have strong odd-even variation due to the blocking mechanism. 
A first order coupling of the neutron pairing tensor to the proton field is not possible in 
a quasiparticle picture. The only possibility is to introduce effective many body forces. 



Zawischa and coworkers [19|], |2(J have introduced a separable four body interaction between 
protons and neutrons within a HFB model and they were able to correctly describe the effect 
for a large range of elements. A four body interaction was indicated by the above mentioned 
requirements for the underlying staggering mechanism, and by the importance of four body 



correlations to the alpha clustering effects [^TJ. Finally, they have shown |p!qj that a 
three body interaction is enough to explain the staggering amplitude. 



In a recent series of papers E3]- [57], an effective four body interaction of monopole 



pairing type, between pairs of protons and neutrons, have been introduced and extensively 



studied. The main consequences of this new interaction, when used in a variational approach 
with BCS wave functions, are (i) the mutual induction of the superfluidity from the neutron 
system to the proton one ||25|| , |[27|| ; (ii) the appearance of local minima in the functional en- 
ergy versus the gap parameters P4"| , [p5f| , which can be interpreted as metastable (superffuid 



isomer) states, possible identified [2~7| with the Oj" state in 152 Sm. The property (i) has been 



successfully used to described the main part of the staggering amplitude for the Pb isotopes 
f27| . This property represents, in fact, the necessarily strong collective coupling between 
the superfluid neutron system and the proton system which become superfluid. A variation 
in the neutron pairing tensor, due to the blocking mechanism in the odd system, induce a 
similar variation in the proton gap. The proton single particle distribution around the Fermi 
level follows this variation and the charge radii of the odd isotopes are lower than the average 
radii of the neighbor even isotopes. This mechanism is in agreement with the qualitative 



explanation given by Zawischa fL8fl . The HFB model include much more matrix elements 
which can connect the long range mean field with the pairing field; this is the reason why 
it includes different type of contributions to the staggering. For example, with three body 
forces only, the main contribution comes from the coupling of the neutron pairing tensor to 
the proton long range mean field [[2tJ which lead to variations of the proton single particle 
wave functions, while with four body forces the main effect is given by the staggering of the 
proton pairing tensor induced by the neutron pairing tensor. 



This debate was a stimulation for the present work. In a previous paper |27]] the Pb 



isotope chain has been studied because for that case no large deformation effects are expected 
due to proton magic number. The staggering amplitude was fairly described for the N ~ 126 
isotopes whose deformations are known to be small [p8[| . A partial agreement was obtained 



for the neutron deficient isotopes possible affected by the deformation |29|| . In this work, 
the Sn isotopes are studied to further establish the contribution of the effective four body 
interaction to the staggering mechanism. The effect of the deformation is also taken into 
account. The Sn isotopes are almost spherical, but their /3 2 's extracted from the BE2 values 
28fl are 3-4 times larger than the similar ones for Pb isotopes. However, their variation 



with the neutron number is rather small and we expect no influence of the shape fluctuation 
to the staggering. It is also interesting to investigate the influence of the particle number 
violation (inherent in any quasiparticle description) to the odd-even staggering effect. To 
this goal we shall compare the results with and without the projection on good number of 
particles. 

In our model, the protons and the neutrons are supposed to move in deformed single 
particle orbits described by a canonical mean field H m ^, and they are affected by a residual 



interaction of the usual pairing type H p<Mr and a new p3| four body interaction of monopole 
pairing type between pairs of protons and neutrons H4 



H= £ {Hf + HD + H A , (1) 

i=p,n 



where 



H? f = E si a t m a w , i = P,n , (2) 

SiCTi 

Hr = -GiP? Pi , P = £a Si _a Si+ , (3) 

Si 

H 4 = -G 4 P;P+P n P p , (4) 

ajf. a is the creation operator of a nucleon in a single particle state with the sign of the angular 
momentum projection on the intrinsic symmetry axis (Xj, Sj representing all the other single 
particle quantum numbers. The reason to preserve only this component from a full four 
body interaction is that all the other components seem to give only a renormalization effect 
to the two body strengths. This conclusion is based on the assumption that the two-body 
pairing strengths are phenomenological parameters, which can incorporate terms depending 
weakly on the extensive properties of the system, like v A (terms like (G p + G Ap xi) still have 
small renormalization effects). As a consequence we keep only that part of the full four-body 
interaction which has a strong dependence on the product of the proton and neutron paring 
tensors. Another reason is its simplicity that permits to obtain complete solutions (in a 
BCS approximation), which exhibit very interesting physical properties. 



The variational solutions of the proposed hamiltonian in the space of BCS wave functions 
are given |25[ by the single particle densities 

<(«l) = ^(l-(+)^=^) , (5) 

the quasiparticle energies 

e Sl = [(E Si -\^ + A^ , (6) 
the constraints for good average number of particles 

Ni = E(l " ~~r~~) (7) 



and the gap equations 



(G P + G,xl) E ^ = 2 (8) 

Sp £ Sp 



(G n + G,xl) E = 2 . (9) 



In the above formulae the E Si are the renormalized single particle energies p5| , Aj is the 
Fermi level for the system i, 

Xi = J2u Si v Si (10) 

Si 

is the pairing tensor for a system with an even number of particles or 

Xi = E u « v si ( n ) 

for a system with an odd number of particles, when the soi level is blocked. The coupling 
constants G p , G n and G4 can eventually be extracted from experimental binding energies 
combination PSJ assuming that they have the following A dependence: 



n — Gp n — Gn n _ C4 . . 

(jr p — — , (jr n — — , <jt 4 — — . : I _ ; 



The gap equations (8) - (9) are strongly coupled by the G4 strength and the collective 
pairing tensor of the complementary system, due to the proposed four body interaction. 
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This coupling has two consequences (i) the proton system, which is in a normal phase for a 
magic number due to the Belyaev condition, can become superfluid due to the contribution 
of the G^xi term; (ii) the staggering of the neutron pairing tensor induced by the blocking 
mechanism is fairly followed by the proton gap A p (see Eq. (8)), the single particle densities 
Vg and finally by the charge radii 

2 

<r2>= ^E <s pl r2 l s p>< • ( 13 ) 

In particular, the charge radii are a little bit larger in an even system, due to the fact that 
a larger proton gap gives larger probabilities for protons to occupy higher single particles 
states with larger single particle radii. 

This simple model has the advantage that it meets the both requirements of Zawischa 
and it can be extended to include quantum fluctuations. It is interesting to try to include 
a 3-body effective force in this type of dynamical calculations (not affecting the details of 
the single particle wave functions but only their superfluid properties). The most simplified 
form of a proton-neutron 3-body effective interaction 

-^3 = ~~G pn 2_j a s p <T p a Sp(T p P n P-n ~~ G np a s n a n a s n cr n Pp Pp ) 

has only the effect of a renormalization of the pairing coupling constants G p (G n ) by the 
quantities G pn N p (G np N n ) proportional with the number of particles. This N n linear de- 
pendence of G p cannot explain any staggering. This represents a clear indication that the 



staggering induced by 3 body forces within the HFB model |22| originates mainly from the 
modification of the shape of the single particle wave functions due to the nonlinear coupling 
of the neutron pairing tensor to the long range mean field. 

One of the drawbacks of the BCS wave functions consists in their nonconservation of the 
particle number. Fluctuations of the particle number give large contributions to quantities 
like the total energy but, they can be important also for such delicate observables like the 
small variation of the charge radii. In order to study such an effect, particle number projected 
BCS wave functions are used 



I* >= P A \BCS >= i- p^'nK + ^< + <-)|- > , (14) 
where N is the particle number operator. The relevant matrix elements can be obtained 

<tf|tf >=R° , (15) 

< ^|a+a s |* >= R\(s) , (16) 



where RJ£ functions can be analytically written, but a recursive relation [ 30 1 



K( Sl ...s m ) = u 2 s R™ +1 ( Sl ...s m s) + v 2 s R™+ 1 1 (s 1 ...s m s) (17) 

is much more suitable for computer calculations. The charge radius is given by 

<r 2 >=-^^< Sp |r 2 | Sp >^ 2 ^( Sp ) . (18) 

The qualitative effect of the particle number projection on the BCS wave functions is that 
the single particle distribution around Fermi surface is squeezed. 

Due to the fact that Sn has a proton magic number and in order to study the effects of 
different parameters, it is interesting to investigate this chain of isotopes in a spherical model. 
For the calculation of the single particle wave functions and energies, a program developed 



by Hird [31] has been used. It has the advantage that the eigenvectors are obtained in terms 
of the tridimensional harmonic oscillator wave function basis for which the matrix elements 
are analytically known. 

The single particle potential is of the Woods Saxon type, allowing /3 2 and (3^ deformations. 
The potential parameters are taken from Ref. page 21 (quoted as Cheprunov in Table 
1 of Ref. For computational reasons (to get double degenerate levels), the calculations 

have been performed with a very small deformation 

/3a = 0.005 , /? 4 = 0.0 . 

The results are summarized in Table 1. The staggering is qualitatively obtained but the 
absolute results deviate from the experimental values. The main difficulty comes from the 
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fact that the slope of variation of the mean square radii versus the mass number A is not 
correctly reproduced. This difficulty has its roots in the shell model results (pairing residual 
interaction neglected) as can be extracted from the third column of Table 1. Possible 
explanations can be: no enough harmonic oscillator shells included, an inaccurate A and 
Z dependence of the Woods Saxon parameters, deformation effects not taken into account. 
The contribution from the higher harmonic oscillator shells has been carefully checked and 
found negligeable. The parametrizations Blomqv - Wahlb, Rost and "universal" from Table 
1 in Ref. |33|] have been checked but no qualitative improvements have been obtained. The 



effect of the deformation will be discussed in the next section. 

The odd-even effect is obtained but the magnitude of the variation and the amplitude of 
the staggering are half the experimental one (see Table 1). In order to obtain this effect a 
fixed set of two and four body pairing strengths has been used 

C p = 33.5 MeV , C n = 20. MeV , C 4 = 2. MeV . 

These parameters can be extracted for every nucleus from the correlation energies Pz, Pn 
and P 4 (defined in Ref. |2"5|) with a relatively complicated procedure described in Ref. P5 



In this work we use the same strengths for the whole chain of isotopes and compare their 
values with the range given by the procedure described in Ref. [^5]. The small slope of the 



radii versus the mass number A, force a relatively large proton pairing strength, C p and a 



small C4 quantity (in Ref. [25[ the experimentally extracted C4 in the rare earth region are 
10 times larger). 



The staggering parameter 10 



2[< r 2 > A+1 -<r 2 > A ] 

7A even = 5 5 > {^) 

< r 2 > A +2 ~ < r 2 > A 

represents an intrinsic measure of the effect. The calculated 7^ (see Table 1) is in very good 
accord with the experiment, indicating that the proposed mechanism is adequate to describe 
the staggering. 

One can try to improve the above results by including the deformation in the model. 
The static deformations for Sn isotopes are unknown but the /?2 values extracted from the 
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B(E2) values can give a fairly good approximation if the deformation is not too small. The 
fa values for the even isotopes can be extracted from the measured B(E2) values according 



to the following formula ||34|| : 



02 = v / < $ > = y/ B ( E2 )^ > (20) 

where Rq is usually taken as 1.4 A 1 / 3 fm. The extracted /?2 values for the Sn isotopes are 
around 0.1 3-4 times larger than the similar values for the Pb isotopes, and their total 
variation is around 0.04 (see Table 2). In order to have an idea of the effect of this variation 
of deformation on the change in the mean square radii one can use the model of a uniformly 
charged deformed nucleus Q 

5 

S < r 2 >= 5 < r 2 > sph + <r 2 > sph —5 < (3 2 > . (21) 

According to this formula a deviation of 0.01 in /? 2 can give, for nuclei with A ~ 120, a 
deviation of 0.01 in 5 < r 2 >. Such a magnitude is just the unit in the staggering. This 
classical result was checked also in a microscopic calculation. 

Another effect of the deformation is that it gives a better slope of variation of the mean 
square radii versus the mass number A. The results are presented in Table 2 and Fig. 1. 
The deformations used in the calculations for the even nuclei have been taken from Ref. 
|28| . For the odd nuclei the interpolated values have been used. Figure 1 indicates that 
the discrepancy in slope has been reduced but the experimental results are not completely 
reproduced. The remaining difference can come from a better A dependence of the Rq 
Woods Saxon parameter. 

The two and four body pairing strengths used in the calculation are 

C p = 30. MeV , C n = 18. MeV , C 4 = 10. MeV . 

The proton pairing strength is not yet decreased very much but, the four body strength 
is near the right magnitude P5fl , The magnitude of the variation of the of the radii 

and the staggering amplitude are now closer to the experimental ones (see Fig. 2). The 
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staggering parameter is smaller than the experimental one (see Fig. 3), indicating that the 
effect is underestimated. 

The effect of particle number projection is shown in Fig. 3 where the staggering param- 
eters with and without projection are compared. The results consist in an increase of ja 
for the projected case in the neutron rich region. This behavior is due to the fact that the 
amplitude of the staggering comes from the amplitude of variation in the proton gap; for the 
neutron rich isotopes the gap is already small in the BCS approximation but, the projection 
shrinks further the proton density around the Fermi level, making the staggering amplitude 
smaller. 

It is interesting to see if the quoted pairing constants are able to reproduce the experi- 
mental pairing energies, 

P z = 1 -[2£{Z - 1, N) - £(Z, N) - £(Z - 2, N)] (22) 

P N = ±[2£(Z,N-l)-£(Z,N)-£(Z,N-2)] , (23) 

where the -£ = B is the binding energy. These quantities are compared to the experiment 
in Fig. 4. The neutron pairing energy is fairly reproduced. The shell effect around A = 
116, connected with a shell closure for N = 64 is also obtained. The experimental values 
for the proton pairing energies, Pz, are larger than 1 MeV. This give some indication that 
the proton system is superfluid to some extent; this effect can be explained by the 4-body 
neutron-proton coupling. The theoretical description of the proton pairing energies is not 
very accurate, especially in the neutron deficient region. This could rise some doubts about 
the agreement obtained for the amplitude of the staggering. The origin of this discrepancy is 
due to a larger C p values necessary to assure proton superfluid properties to all even isotopes 
in the chain. However, if one wants to describe only the neutron deficient part of the isotope 
chain with a fixed set of strengths, one could do this with a lower value of C p and, as 
a consequence, with a better description of the proton pairing energies. A more accurate 
description of the neutron number dependence for the slope of the radii and for the G p , G n , 
G 4 strengths (see Eq. (12)) could solve this discrepancy. 
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Finally, it is interesting to compare the absolute rms radii with the values experimentally 
extracted || . The results are given in the last column of Table 2. They are in good agreement 
with the experiment for the neutron deficient isotopes (compare to the data in the prelast 
column of the table). Similar discrepancy have been obtained in the more sophisticated HF 
calculations for Te isotopes (see Table XIV in Ref. [|J). 

In conclusion, the staggering of the nuclear charge radii of the Sn isotopes has been in- 
vestigated in a simple superfluid model with an effective four body interaction of monopole 
pairing type included. The intrinsic effect has been obtained in a simple quasispherical ap- 
proximation but, the absolute results deviate from the experimental values due to the weaker 
dependence of the slope of the charge radii versus the neutron number. The contribution of 
the particle number projection and of the deformation have been studied, their effects lead- 
ing to better absolute results. This simple model has two interesting features: (i) it fulfill the 
qualitative physical requirements of the staggering mechanism discussed by Zawischa; (ii) 
it can be extended to include quantum fluctuations like particle number projection, RPA 
pairing and quadrupole vibrations in the ground state or fluctuations in the gauge space 
connected with the gaps. Its parameters have known physical meaning and its ability to 
describe the staggering data, make it a good candidate model to check the contributions of 
different mechanisms to the odd-even staggering of the charge radii. 
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Table Captions 



Table 1 Relative and the variation of the mean squared charge radii (in fm 2 ) for Sn iso- 
topes. Experimantal values are taken from Ref . . Theoretical values obtained in the 
quasispherical model. Last two columns present results for the staggering parameter, 
Eq. ©. 



Table 2 Charge radii results with deformations included. First three columns give the 
relative mean square charge radii in fm 2 . /?2 values are taken from Ref. [28|. Last two 
columns show the absolute charge radii in fm; experimental values from Ref. M. 
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Figure Captions 

Figure 1 Relative mean squared charge radii for Sn isotopes. Experimental values (squares 
and solid line) taken from Ref. 0. Theoretical values (circles and solid lines) are cal- 
culated with four body forces and deformation included. Theoretical values calculated 
with projected wave functions are denoted by pluses and are connected by dashed 
lines. All lines are drawn to guide the eye. 

Figure 2 Staggering of the mean squared charge radii. The significance of signs and lines 
is the same as in the caption to Fig. 1. 

Figure 3 Theoretical staggering parameter compared with the experimental one. The 
significance of signs and lines is the same as in the caption to Fig. 1. 

Figure 4 Neutrons (left) and protons (right) pairing energies given by Eqs. (|23|) and (p2|). 
Experimental values marked by squares and solid line. Theoretical values (circles and 
dashed lines) are calculated with four body forces and deformation included. All lines 
are drawn to guide the eye. 
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A 


Exp. 


Shell model 


4N model 


5 < r 2 >a,a-i 




{1a)an 


110 


-0.638 


-0.439 


-0.360 




0.736 


0.876 


111 


-0.586 


-0.388 


-0.323 


0.037 






112 


-0.497 


-0.339 


-0.275 


0.048 


0.908 


0.806 


113 


-0.438 


-0.293 


-0.244 


0.031 






114 


-0.367 


-0.246 


-0.198 


0.046 


0.687 


0.706 


115 


-0.322 


-0.202 


-0.173 


0.025 






116 


-0.236 


-0.159 


-0.127 


0.046 


0.758 


0.718 


117 


-0.189 


-0.118 


-0.103 


0.024 






118 


-0.112 


-0.077 


-0.060 


0.043 


0.750 


0.707 


119 


-0.070 


-0.038 


-0.039 


0.021 






120 


0.0 


0.0 


0.0 


0.039 


0.832 


0.732 


121 


0.042 


0.038 


0.019 


0.019 






122 


0.101 


0.074 


0.052 


0.033 


0.725 


0.723 


123 


0.134 


0.109 


0.069 


0.017 






124 


0.192 


0.144 


0.099 


0.030 






125 


0.225 


0.178 


0.112 


0.013 
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A 


Exp. 


4N model 


Projected 


fa 


rms exp. 


rms 4N model 


110 


-0.638 


-0.504 


-0.512 


0.126 




4.5820 


111 


-0.586 


-0.462 


-0.467 


0.1245 




4.5866 


112 


-0.497 


-0.386 


-0.390 


0.1227 


4.5958 


4.5949 


113 


-0.438 


-0.329 


-0.328 


0.1208 




4.6011 


114 


-0.367 
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